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1 Introduction 

Astrophysical data show that the equation of state parameter lo lies near a; = — 1 
and possibly uj < —1, leading to an accelerated universe Some present data 
seem to favor an evolving dark energy with lo less than —1 at present epoch from 
u; > — 1 in the near past [2]. So we may assume that the universe is dominated by 
a perfect fluid for which w < — 1, dubbed phantom energy mu. This description 
of the universe may contain finite-time future singularity, accompanied with 
divergence of dark energy density, called Big Rip P). Depending on the phantom 
potential, this model may also lead to universes expanding for ever, e.g., universe 
tending to a de Sitter space time |^. The effects of gravitational back-reactions 
can also counteract that of phantom energy and can become large enough to 
terminate the phantom dominated phase before the big rip |^. 

Thermodynamics of the expanding universe has also been the subject of sev¬ 
eral studies Phantom thermodynamics looks 

leading to negative entropy of the universe m or to appearance of negative tem¬ 
peratures HU. In accelerated expanding universe, besides the normal entropy, 
a cosmological horizon entropy can also be considered. One can investigate the 
conditions for which the generalized second law of thermodynamics (GSL) holds 
[naHH. In these cases GSL asserts that the sum of the horizon entropy, and 
the normal entropy of the fluid is an increasing function of time. In m the 
change in event-horizon area in cosmological models that depart slightly from 
de Sitter space was investigated, and it was shown that the area and conse¬ 
quently the (de Sitter) horizon entropy are non decreasing functions of time. In 
the presence of a viscous fluid, there was found that GSL was satisfied provided 
that the temperature of the fluid was equal to or lower than de Sitter horizon 
temperature, Tds = H/{2Tr), where H is the Hubble parameter. In HSl the GSL 
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for two specific examples of phantom dominated universe was discussed. In the 
first example it was shown that for a time independent parameter of state, the 
total entropy is a constant. In the second example using a power law potential 
and the slow climb approximation, it was shown that GSL was satisfied. In the 
absence of a well defined Hawking temperature for cosmological horizon, the 
temperature of the phantom fluid in m was assumed to be the same as the de 
Sitter temperature. 

The phantom universe, has a future event horizon (cosmological horizon), 
Rfi- In this model the Hubble parameter is an increasing function of time so 
we expect that the horizon entropy, in contrast to the model discussed in nm, 
decreases with time. 

In this paper we consider the phantom dominated era, and prove that the 
future event horizon area is a non-increasing function of time. It is shown that 
one necessary condition for satisfying GSL is the positivity of the temperature. 
If the temperature is assumed to be proportional to the de Sitter temperature 


m 


T = 




( 1 ) 


an inequality in terms of Hubble parameter, iJ, the future event horizon, Rh, and 
the parameter b, will be obtained. By solving this inequality one can determine 
the dynamics of Rh, H, and also the range of b. We show that GSL is satisfied 
for a wide range of models, provided that 0 < b < 1. We also study the condition 
of validity of GSL in the transition from quintessence regime to the phantom 
dominated era. At the end, two specific examples are used to emphasize our 
general results. 

We use the units h=c = G = kB = l- 


2 Phantom thermodynamics 


The spatially flat FRW universe in comoving coordinates {t, x, y, z), is described 
by the metric 

ds^ = —dt^ + a^(t){dx‘^ + dy'^ dz"^). (2) 

The relative expansion velocity of the universe, in terms of the scale factor a(t), 
is given by the Hubble parameter H = d/a. The over dot indicates derivative 
with respect to the comoving time t. We assume that the universe is filled 
with a perfect fluid with the energy density p, and the pressure P. The fluid 
state parameter uj is defined through the equation of state P = ujp. Einstein 
equations are 


dH 

dt 

= 


—47r(P -I- p) 
8tt 


The above equations result 


UJ = —1 — 


2H 


(3) 

(4) 


For an accelerating universe we have H + > 0, which, in terms of uj, can be 

expressed as a; < — 4. The corresponding cosmological fluid is then consisted of 
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some sort of energy density which has a negative pressure known as dark energy. 
Matter with a; < — 1 is dubbed phantom energy. By considering 0 , one can 
see that for phantom fluid: H > 0, which using © results P + p < 0. 

There are various approaches to study the dark energy. One approach is 
introducing scalar helds. The behavior of io in the phantom regime, uj < —1, 
can be related to the presence of a phantom scalar held, 4>, with a wrong sign 
kinetic term |S]. Depending on the form of the phantom scalar held potential, 
different solutions such as asymptotic de Sitter, big rip, and so on may be 
obtained p]. To describe ui > —1, or quintessence regime, a normal scalar held, 
(T, known as quintessence scalar held, can be used cni. The phase transitions 
in inhationary models can be investigated using hybrid models. These models, 
at least, are composed of two scalar helds EOI, e.g., one of these models is the 
quintom model which in order to describe the transition from a;>—ltoaj<—1 
regime, assumes that the cosmological huid, besides the ordinary matter and 
radiation, is consisted of a quintessence and a phantom scalar held EH 

Another method to study the acceleration of the universe is to use a running 
cosmological constant based on the principles of quantum held theory (spe¬ 
cially on the renormalization group) which can mimic the behavior expected 
for quintessence and phantom -like representations of the dark energy. In this 
method despite of the absence of scalar helds, an effective equation state like 
0 can be obtained m 

In this paper we will use the equations 0 , 0 , and assume that the universe 
is hlled with a perfect huid in phantom phase in the sense that to < —1, but our 
results will be independent of the scalar helds or other origins of the phantom 
energy. 

Consider the era when phantom phase is dominated i.e, H > 0. We assume 
that a{t) —>■ oo when t ts, so the scale factor will diverge for a future value 
on the world time: ts- The radius of the observer’s (future) event horizon is 


Rh 




Using 0one can see that Rh satishes the following equations 


(5) 


Rh = HRh — 1 , 

iih = {H + H^)Rh-H. (6) 


For de Sitter space time, tg = oo, Rh = ^/H, and the eq.0 continues to hold. 

Theorem: Rh is a non-increasing function of time, Rh < 0. 

Proof: We have H > 0, therefore d > of/a. This inequality can be integrated 
from t to ts 




(7) 


resulting 


I 

a{t) 


1 

a{ts) 



( 8 ) 


Suppose that is bounded from below by a positive number e. In our 

(super)accelerating model of universe, d > 0 yields 


I 1 

a{t) a{ts) 


t < ts. 


(9) 
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Therefore 


( 10 ) 



But (lion conflicts with the assumption of having a future horizon©. Therefore 
l/a{ts) = 0, and using ©, we obtain HR^ < 1. Using this result and © , we 
arrive at < 0. 

The future event horizon area is We assume that this horizon area 

represents a true contribution to the total entropy. The corresponding entropy 
is proposed to be M 

Sh = ttRI ( 11 ) 

In this way Sh = 2 ttR hRu which is negative and shows that Sh decreases with 
time. This is in contrast to the assumption of m, upon which it was argued 
that the holographic dark energy has no phantom-like behavior. 

The entropy of the phantom fluid inside the cosmological horizon of a co¬ 
moving observer is related to the energy and the pressure via the first law of 
thermodynamics 


TdS = dE + PdV = {P + p)dV + Vdp. (12) 

In terms of H this law can be written as 

TdS = -^^dV + Vdp. (13) 

Air at 

Using V = ^TrRh and Pl^ = (87r/3)p, we obtain 

TdS = -HRldRh + HRldH. (14) 

Then using , we arrive at 

TS = HRl- (15) 

For phantom fluid H > Q, S > Q, provided we take T > 0. GSL states that the 
sum of the ordinary entropy and the horizon entropy cannot decrease with time 

S + Sh>Q. (16) 

In our model this law yields 


-Rh + 2^Rh > 0. (17) 

Note that in phantom regime, H > 0 and Rh < 0, so, a necessary condition 
that the GSL to be satisfied is positivity of the phantom fluid temperature 
T > 0. Therefore the GSL is violated in phantom models with positive entropy 
and negative temperature CHI For a discussion about this subject see El, m 
The GSL asserts that the entropy of an isolated system cannot decrease, but 
does not determine the sign of the entropy. But if we assume that at the finite 
time t = ts(i.e. when Rh — 0) we have S = 0 d, and the GSL is assumed 
to be satisfied in an era before ts, then the entropy will be negative in this 
epoch, which agrees to d and El in which a negative entropy for a specific 
model of phantom thermodynamics with positive temperature has been derived. 
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However in the above mentioned papers i.e., ini, m and ng, the role of the 
future event horizon has been ignored. 

If the temperature is taken to be T = |f-, GSL results 


b 


< 

- Hr^- 


(18) 


For a de Sitter space-time Rh = jj, which results & < 1, so we expect that 
in a phantom model which is small perturbed around de sitter space the fluid 
temperature is equal to or less than the de Sitter temperature. 

If the total entropy is increasing in the universe expansion, 


S + Sh > 0, 

equation results 

{mRh) > 0 . 

For 6 = 1, we must have {HRh) > 0, or using © 

iih > 0 . 


(19) 

( 20 ) 

( 21 ) 


Here Rt is an increasing negative function of time. When Rh 0“, Rh tends 
to the bubble radius 1/H (see equation 0). 

Now let us assume that the expansion of the universe is a reversible adiabatic 
process such that 

S + Sh=0. ( 22 ) 

From we obtain HR’^ = A where A is a constant. 

For 6=1, using this results Rh = X — 1 and using the fact that Rh{ts) = 0 , 
for A ^ 1, we obtain 

Rh = {l-X){ts-t). (23) 

Note that Rh <0, therefore A < 1. w is given by 


1 2 

uj = -< — 1 . 

3 3A 


(24) 


If we take tg = oo, Rh becomes infinite as claimed in |g. If A = 1, the space-time 
is de Sitter, Rh — ^ and w = — 1. 

For 6^1, (ISg and the first equation of lead to the following equation 


Rh - XRl-’’ + 1 = 0, 


(25) 


with solutions satisfying 

i?/.d>(A^^^^ 1, ^) 


(26) 


$, is the Lerchiphi function. The condition Rh < 0 yields A < but 

Rh{ts) = 0, so to have A 7 ^ 0, 6 must satisfy the inequality 6 < 1, i.e, T < 

Rh = 0 occurs at Rh = A'^^, which is a branch point of $. At t = tg we have 
= 0 as expected. See fig. O- 

The above results are general and do not depend on the previous phase of 
the phantom universe. As it was mentioned in the introduction, based on recent 
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Figure 1: ts — t a.s a. function of Rh, for b = 0.9, A = 1 (continuous line) and 
A = 1.1 (points). 


astrophysical data, one can consider the possibility of the transition from uj > —1 
to w < — 1. So it may be interesting to see how our results may be restricted, 
if the GSL still hold in the transition epoch. By transition epoch we mean a 
little neighborhood around the transition time: t*. Like the phantom universe, 
for a quintessence universe we can consider the cosmological event horizon via 
©• The relation sill valid. By considering that H < 0, we deduce that a 
necessary condition for satisfying GSL in such universe is Rh > 0. But in the 
phantom universe, Rh is non-increasing. Therefore for a continuous Rh (which 
may be the result of continuity of H and Rh, see (0), we must have Rh{tt) = 0. 
Using ®, this results HRh = 1 at it- At the transition time, we have also 
H{tt) = 0. In this way using lUTll and dia, we obtain S{tt) = Shitt) = 0, 
therefore the total entropy is differentiable and thereby continuous at tj. For 
t < Sh is an increasing (decreasing) function, while the fluid entropy is a 

decreasing (increasing) function of time. If the entropy is taken to be negative 
in the beginning of phantom dominated era (which is not a consequence of GSL 
which considers only the time derivative of the entropy) and the GSL is assumed 
to be valid in the transition epoch, then we must have a negative entropy in the 
quintessence era too. 

3 Examples 

To illustrate our general results, let us consider two examples. As a first example 
consider the phantom universe of pole-like type 

a(t) = ao(G - t)"”, (27) 

t < ts and n > 0 is a positive real number, oq is a real positive constant. 
H > 0, and uj = —1 —2/(3n) < —1 is a time independent constant. In this model 
HRh = n/{l+n), and Rh = —l/(n-|-l) < 0. This is consistent with the theorem 
imposed after the eq. ®. Following the discussion after the eq. ( 1 ^ . and by 
taking A = n/(n + l), this model corresponds to a reversible adiabatic expansion 
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of an universe with the phantom temperature T = H/ (27r) in agreement with 
the results of m 

To show the role of the parameter b in validity of GSL, as an another example, 
consider the model m 

“ = (28) 
where oq and n are two positive real constants. In this model 


H 

H 


ts 

""tits - t) 
nts{2t - ts) 
t^ts-t)^ ■ 


(29) 


The condition 77 > 0 is satisfied when ts < 2t. lo is time dependent 

2{2t-ts) 

UJ — — 1 — ---. 

3nts 

When 2t > ts, i.e., in phantom dominated era, we have uj < —1. 

In this model the future event horizon is 

Rh = ts{x — l)~^ J {u — 

where x := ts/t, \ < x <2. We have 

Rfi =—1 + nx^{x — J {u—l)^u~^du. 


(30) 


(31) 


(32) 


Note that from the above equation one can show that R^ < 0, in agreement 
with our previous results. This can be verified by considering that for x <2 we 
have 


n 




l)-((n-l)u + 2) ^^^ 


{x- l)”+i 


(33) 


Which can be verified by noting that the integrands satisfy the above inequality 
for all u belonging to [1, 2]. 

If GSL is respected, we must have (H^Rh) > 0, where c := 1/6. This 
condition can be written as 

(x — l)“"“^((n — c)x^ + 2cx) J {u — l)^u~‘^du > 1. (34) 


For 1 < X < 2, n > 0, and c > 0 (which follows from the discussion after the 
eq. 113), we have (n — c)x^ + 2cx > 0, therefore, (IHUl reduces to 



l)"u ^du > 


(x- l)”+i 
(n — c)x^ + 2cx 


(35) 


The right side of the above inequality can be written as a definite integral 

(1 — n)(c — n)v? + 2(n + cn — c)u + 2c 


/ 


{u—l)^u ^du > I (u—l)"u 


((c — n)u — 2c) ^ 


du, 

(36) 
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which can be rewritten as 



(m — 1)”m ^ 
((c — n)u — 2c) 


qdu > 0. 


(37) 


q is defined through 


q := Au^ + 2Bu + D, 


(38) 


A := (c — l)(c — n), B := (n{c — 1) + c — 2c^) and D := 4c^ — 2c. If for all u 
in the interval u S [l,a;] C [ 1 , 2 ], q > (<) 0 , then the inequality holds (does 
not hold). At u = 1, q{u = 1) = (c — l)(n + c) and at u = 2, q is negative: 
q{u = 2) = —2c. Therefore a necessary condition for validity of GSL is that q 
must have two roots such that at least one of them lies in [1,2]. Otherwise q 
will be negative in this interval. 

If n and c satisfy 


> 


c^(l — 2c) 

(c-l)2 


(39) 


q has two roots. 

For A > 0, g has a root l<r<2atr = {—B — y/ B^ — AD)/A, provided 
that c > 1. The other root takes place at r' = {—B + y/B"^ — AD)/A > 2. 

For A < 0, g has a root at 1 < r = {—B — y/B^ — AD)/A < 2, provided that 
c > 1. The other root takes place at r' = {—B + y/B‘^ — AD)/A < 1. 

In these cases we have g(u = 1) > 0, hence g is positive for all values of u 
satisfying 1 < u < r and negative for all values of u in the interval r < u < 2. 
For c < 1, neither of the roots of g lies in the interval[l, 2] and the sign of g is 
negative. 

For A = 0 (in this case g(u) is a line), we have either c = 1 or c = n. 
For c = 1, q{u = 1) = 0 and g is negative in (1,2]. When c = n, we have 
q{u = 1 ) = 2 c(c— 1 ) and q{u = 2 ) = — 2 c, so if c < 1 , g will be negative in [ 1 , 2 ]. 

Using these results we conclude that g has a root in [1,2], provided that 
c > 1. Other values of c lead to negative values of g for all u G [1,2]. Therefore 
T < H/{2tt) is a necessary condition for validity of GSL. 

Despite the above result, c > 1 is not an enough condition for validity of 
GSL, specially in the early period in which the phantom became dominated 
as is illustrated in fig. ©■ Using our arguments we can only claim that in the 
period 1 < a: < r, (in an interval before the big rip)the GSL holds. But note 
that for large c or large n we have r ~ 2, and GSL is applicable approximately 
in the whole region of the interval [ 1 , 2 ]. 

The transition from H < 0 to H >0, occurs at x = 2, or t = where H = 
0, which results 5(U) = 0. As we have discussed in the last part of the previous 
section, a necessary condition for satisfying GSL in t is Rh{x = 2) = 0 

which using 123), implies 


4n y {u— l)"u '^du = 1. (40) 


The above integral can be written in terms of digamma functions, resulting 


2 (n + l) 


1 1 , T 

—-- + - n 4- 

2 (n+l) 2 



An 


1. (41) 
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Figure 2: /“(«—!)”« — (a; — l)"+^((n — c)a;^ + 2ca;) as a function of a; 

for n = 1 and c = 2. 

Following our previous discussion, we expect that GSL remains valid in the 
region near a; = 2 for large n. By asymptotic expanding of digamma functions, 
the above equation, becomes 



(42) 


which is valid for large n. Therefore for large n, the GSL may be true in the 
transition epoch. The condition 123, in a short time before the transition: 
t = tt — e = ts/2 — e, becomes 



(43) 


which based on our previous discussion about the enough condition for validity 
of GSL, is correct for large n. Hence the GSL is respected in the transition 
epoch only for large values of n. 

4 Summary 

In this paper we discussed the constraints and conditions imposed on cosmo¬ 
logical future horizon Rh, Hubble parameter H, and the temperature T, in a 
phantom dominated universe in order to satisfy the GSL. We proved that the 
future event horizon is a non-increasing function(see the theorem after eq.®) 
and thereby the corresponding horizon entropy is also a non-increasing func¬ 
tion of comoving time. To obtain this result we assumed that one can relate 
an entropy to the future cosmological horizon, in the same way as black-hole 
event horizons, despite their different structures cnmiiiini. Using the assump¬ 
tion that the (super) accelerated universe is filled with a perfect fluid satisfying 
the first law of thermodynamics (eq. HI 211 1. we obtained the total entropy and 
found that the positivity of the temperature is a necessary condition for GSL 
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(see eq. C3). We obtained also an inequality in terms of Rh , H^ and T (see 
eqs. (HU), (HOI) and HU). This inequality becomes an equation for the adiabatic 
reversible expansion of the (super) accelerated universe. We solved this equa¬ 
tion for Rh, for de Sitter and non-de Sitter temperatures (see eqs. (HU and Hnj). 
We showed that in various cases T must be less than or equal to de Sitter tem¬ 
perature. By studying the influence of the transition from the quintessence to 
phantom dominated universe on the GSL , we deduced that the time derivative 
of the future event horizon and the entropy must be zero at the transition time. 
However if we adopt that the GSL is valid in the transition epoch and also as¬ 
sume that the entropy is negative in the beginning of the phantom era (which is 
not a consequence of GSL), a negative entropy in the quintessence regime will 
be deduced, which is not expected. 
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